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Abstract 

We present a calculation of the i3°i? -mixing including Next-to-Leading Order (NLO) QCD cor- 
rections within the Two Higgs Doublet Model (2HDM). The QCD corrections at NLO are contained 
in the factor denoted by 772 which modifies the result obtained at the lowest order of perturbation 
theory. In the Standard Model case, we confirm the results for 772 obtained by Buras, Jamin and 
Weisz [U . The factor 772 is gauge and renormalization prescription invariant and it does not depend 
on the infrared behaviour of the theory, which constitutes an important test of the calculations. 
The NLO-calculations within the 2HDM enhance the LO-result up to 18%, which affects the 
correlation between Mh and Vtd- 

PACS-numbers: 12.15.Ff; 12.60.Fr; 14.80.Cp; 12.38.Bx 

Key-words: Mixing; Oscillations; 2HDM; CP-violation; CKM-matrix 



1 Introduction 

The 5-meson plays an important role in present day physics. i?°I? -mixing as well as 5-meson 
decays can be used to determine the CKM-elements and to investigate CP- violation within the 
Standard Model (SM). Two large experiments, namely the BABAR-Collaboration at SLAC 
and the Belle Collaboration at KEK |J, will start taking data in the near future. The major 
subject of research in those experiments is the determination of the CKM-matrix elements. 
This may also provide first hints on physics beyond the Standard Model. 

In general terms 5°B°-mixing is a Flavour Changing Neutral Current (FCNC) process gen- 
erated through weak interactions. In the SM, the process is generated at the lowest order of 
perturbation theory via the box diagrams displayed in Fig. [1|. Due to the large difference in 
the masses of heavy particles in the box (W ± , H ± and t-quarks) and much lighter external 
particles {b and d-quarks), it is possible to disentangle long and short distance effects. This can 
be accomplished within an effective theory. The proper separation of short and long distance 
QCD-effects at Next-to-Leading Order was first presented for the SM in We are utilizing 
here the same framework of renormalization group improved perturbation theory. 



Our purpose is to include two Higgs-boson doublets as an extended Higgs-boson sector. Such 
models are usually called Two Higgs Doublet Models (2HDM) 0, |5]. Two particular models 
which are usually considered (so-called Model I and Model II) differ in the couplings of the 
charged Higgs-bosons to fermions (see for example [[|): 

^Hud = 9 % Ud [m u cot{(3)P L -m d cot(f3)P R ] , (1) 

£>H*d = ^^[m u cot(j3)P L + m d t&n{P)P R } . (2) 
y2m w 

Here gw is the weak coupling constant, m w is the VT-mass, and the projectors Pl,r are defined 
by 

ft*=ip. (3) 

The subscripts u and d denote up- and down-type quarks, respectively. 

As we neglect all quark masses except m t , both models are identical in our case. Thus our 
calculations are valid for all values of tan j3 for the Model I, but only for tan /3 << mt /mj, in the 
Model II. Looking at the Feynman rules for the Higgs-quark-quark coupling, it is easy to verify 
that in Model II, the term proportional to m& becomes important for large values of tan/3. Due 
to the projectors in the vertex, we would end up with different operators for the tan/? and cot /3 
parts. However, it should be stressed here that Model II is more interesting, because this is the 
choice favoured by the Minimal Supersymmetric Standard Model ||, |5j . Large values of tan (3 
will be considered in a future publication. 

We continue our discussion with a short review of the Leading Order (LO) calculations. The 
basic idea is rather simple. First, one has to calculate the box diagrams in Fig .[I]. They can 
be evaluated in the Feynman-'t Hooft gauge for the iy-boson. This leaves us with a physical 
W-boson and a would-be Goldstone boson (unphysical scalar Higgs boson) of the same mass 
and with couplings proportional to the quark masses. As all quark masses except the top-quark 
mass are equal to zero to a good approximation, t is the only active flavour in boxes containing 
Higgs particles or unphysical scalars. In the remaining box diagrams, the effects of u- and c- 
quarks are taken into account by the GIM mechanism. Evaluating the box diagrams in this 
framework leads to the well known Inami-Lim functions 0-[lC . 



The scaling behaviour from the matching scale /j,q = mw of the full and effective theory 
down to a lower scale is then determined by the one-loop QCD correction to the effective 
vertex generated by the previous procedure of integrating out the internal heavy degrees of 
freedom. The renormalization group equation gives us the factor rj LO describing the effect of 
scaling [|l]J . 

The effective Hamiltonian for 5°I? -mixing at LO reads 



1 G 2 

H eS = -. — f m 2 w (V td V t * b ) 2 f] LO S(x w ,x H ) O ll , (4) 



4 7T 2 

where 



6 LL = d 1 ,P L b d>fP L b (5) 
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and 



S(x w ,x H ) = S ww (x w ) + 2 S wh (x w ,xh) + S H h(xh) ■ 



(6) 



The LO Inami-Lim functions Sww( x w)> Swh(xw, %h) and Shh(xh) can be found in Appendix 
A. The arguments of the Inami-Lim-f unctions will be denoted xw,h = m t/ m wH- 

The factor r\ LO is determined by the /3-function |12] and the anomalous dimension 7 to be 
obtained by evaluating the diagrams of Fig. 0. It reads 



" a s (m w ) ' 


7 (0)/(2/3 ) 


' a s (m w ) ' 


6/23 


. a s (m B ) _ 




. ot s {m B ) _ 





Vlo 



For later convenience we define t]lo = otsiynw 



0.85 . 



(7) 



i (6/23) 



To evaluate the matrix element of H c g we have to employ 

(B°\6 LL \B°) = I B B (//) f B m% , 

where f B is the 5-meson decay constant and B B ([i) parameterizes deviations from the vacuum 
insertion approximation. 



(8) 



2 Explicit QCD-Corrections 

The perturbative result for 0(as) corrections in 2HDM will be presented in this section. The 
results can be obtained by evaluation of the diagrams of Figs. |3] and ||. We performed our 
calculation in an arbitrary covariant £-gauge for the gluon and we employed the Feynman-'t 
Hooft gauge for the ty-boson. 

As one easily notices, the diagrams a, b and f - i have the "octet" -structure T a ®T a , whereas 
the diagrams c - e and j have "singlet" -structure 1 ® 1 in colour space. The double penguin 
diagram k contributes to the considered process at the order of 0(1/ m^) which is negligible 
for our purpose. 

Diagrams containing vertex- and self-energy corrections (diagrams c, d and e) lead to UV- 
divergent integrals and hence they have to be regularized and renormalized. We are using 
dimensional regularization with anticommuting 75. This corresponds to the NDR-scheme 



"naive" dimensional regularization scheme) [O, 14]. We performed the necessary renormaliza- 



tion in the MS'-scheme |T3| , |T5| . All other diagrams are UV-finite. Furthermore, one notices 
that the diagrams g - j contain infrared divergences. To deal with them we keep the external 
quark masses whenever necessary. As we will see, the external quark masses do not affect the 
final result for the Wilson-coefficient. This justifies the method we have chosen. 

We could also treat the infrared divergences with the method of dimensional regularization 
16| , but our intention was to verify the results provided by Buras and collaborators in |1[] and 
to extend these calculations to the 2HDM. 

The (9(as)-corrections to the Hamiltonian (^) have the following structure 

(Atf eff > = - m 2 w ( Vtd V; b f U(x w , x H ) , (9) 



3 



where 



U(x w , x H )=J2 (c F i®i <$\x w , x H ) +f a ®f a $\x w , x H ) ) (O k ) 



a AS), 



(10) 



with k E {LL, 1, 2, 3}. C F is the colour-factor defined by C F = (N* - 1)/{2N C ) and N c is the 
number of colours. 
The operators Oi read 



d 3 



d P L b] [d P L b] - j [d P L b] [d P L b] + (P L 



Pr) 



dP L b 



dP R b 



+ 



dPnb 



dP L b 



d lfl P L b d-fP R b +(Pl^Pr) ■ 



(11) 
(12) 

(13) 
(14) 



The operator 0\ steems from diagrams g and h in Fig. f|, whereas the operators O2 and O3 
follow from diagrams i and j, respectively. Since the relevant operator Oll is self-conjugate 
under Fierz-transformations, we will use 



f a ®f a = c A i®i = ^(\-^\\, 1 



(15) 



to retain only one operator. Nevertheless, to keep the calculations transparent, we will abandon 
the distinction between octet and singlet at the very end only. 
The coefficient functions in equation fllOD can be decomposed as 



id) 



Xj M \x w ) + Xj""'{xw,x H ) ■ 



(16) 



The functions xf' S i x w) were already given in ref. M. We have recalculated them and we 



confirm these Standard Model results 

xif M) M = L^ SM \x w ) + 



2£ +2£g m +2£\n(x tl0 ) + 6\n(x tl0 ) x w 



d 



dx 



w 



X { ll M \xw) = L^ SM \x w ) + [2i + (3 + Hx b x d ) + 2£g IR ] S ww (x w ) 



X 



(8,SM) 



{x w ) 



-(3 + £) S ww (x w ) 



xf SM) M 



; -2 X ^ SM \x w ) =-(3 + 
with all other x equal to zero. 



m b m d (x d \ 

In I — ) S ww {x w ) 

\ Xfo , 



ml -ml 



Sww{%w) , 

(17) 
(18) 
(19) 

(20) 
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We have introduced the following abbreviations: 



X W,H 

9ir 



m 



in 



%d,b 



d,b 



x 



/'o 



W,H 



III 



/<() 



W 



III 



2 ' 
W 



x d \n(x d ) - x b \n(x b ) 



(21) 
(22) 



x d - x b 



Here, mt stands for the top quark mass renormalized at the scale fio in the MS-scheme. 
The function Xll^ reads 



Xll\ x w,x h ) = L {1 ' H \x w ,x H ) 

2£ +2tg IR + 2£m(:r Mo ) + 6 1m> Mo 

•(2 Swh(xw, %h) + Shh{xw-, x h)) ■ 



v _d_ 

i=H,W OXi 



(23) 



The analytical expressions for the functions Sww, S\vh, Shh and L^' SM,H ' are listed in Ap- 
pendix A. Note that powers of tan (3 enter in the definition of the various 5" and L indicating 
the number of Higgs-bosons involved. 

The remaining newly calculated functions Xll ^ > 
spectively from x^ @, X?' SM) (0), xf SM) © and X ? ' 5M) © by changing L^^(^) 



xi 8 '^\ X2 8 '^^ and xls* can be obtained re- 



to L (i ' h) (xi4/,xh) and to (2 S WH (x w ,x H ) + S HH {x H )). 

We have obtained the function L^' SM ' H ^ by using computer algebra systems, especially the 
program FORM to evaluate the Dirac-structure and Mathematica 3.0 |T7[ for summarizing all 
terms. Maple 



lq| was used for some integrations. 



Note that the results of the particular diagrams exhibit non-trivial dependence on the gauge- 
parameter £ and the IR-masses, but the functions L^' SM ' HS) are gauge as well as renormalization 
prescription invariant. They do not depend on the external quark masses either. Terms de- 
pending on £, x b , x d , or x M0 are proportional to the LO-Inami-Lim-functions Sww, Swh or Shh- 
This reflects the exact factorization of long and short distance effects. 



3 Matching and running 

We start this section with the observation that Eq. (|I0"D is obviously unphysical because the 
coefficient functions Xk are gauge dependent. Actually this is nothing but an artifact of 
the specific way we have regularized the IR-divergences. Other regularization schemes being 
not based on the use of small masses for the external quarks would change the obtained result. 
Dimensional regularization or a gluon mass, for example, would leave us with only one operator 
Oll at this stage. To take proper account of this fact, we have to evaluate the matrix element 
of the physical operator Oll u P to order 0(a>s) using the same IR-regularization prescriptions 
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as before. Yet the one-loop amplitude of Oll as given by Fig. |2| results in the same unphysical 
operators Oi 2 3 with the same coefficients. 



fr(ft))lloop = (Octree + E [^A^M + CUxg^)] 1 ® i<4>tree 



4tt 



where the sum over k = LL, 1, 2, 3 is understood. The functions XAkiVo) are: 



(24) 



Xa?ll = -3 + 2fln(x M ) + 2f-2f 



[x d ln(x d ) -x b ln(x b )} 



(25) 



(8) 

Xa,ll 



v (l) 

AA,3 
(8) 

Xa.i 



6 ln(x M ) - 5 + 2£ + (3 + £) ln(x d x b ) - 2f - 



1 (8) 

-(3 + 



(3 + m d m b /x d 



2 — 



In 



x b 



[x d ln(x d ) - x b ln(x b )} , (26) 

(27) 
(28) 



Eq. (|2J) allows us to write for the matrix element of the effective Hamiltonian up to order 
0(as) in the following form: 



1G 2 

4 7T 2 



r(l)\ 



~a ~~~f~ m w (VtdVtl) 2 C LL {fi ) (O LL (n )) Uoap 



with (O LL (yU ))iioo P as given in Eq. (|3p and C LL (/x ) reads: 



Cll(ho) = S(x w ,x H ) + 



47T 



D(x w ,xh,x, 



IM) i 



(29) 



(30) 



The change due to the inclusion of the charged physical Higgs results in replacing Sww{ x w) 
and D SM {x w ) by 



) = D SM (x w ,x tM) ) + D H (x w ,x H ,Xw) ' 

where 

d 



(31) 
(32) 



D SM (x w ,xJ = C F \L^ SM \x w )+ ^61n(x 



<9x 



ir 



+ 3 



61n(x w ) + 5 



(33) 



6 



D H (x H ) = C F {l {1 > h \x w ,x h ) 



+ 



61n(^ ) J2 



-H,W 



d_ 

dxi 



+ 3 



2 Swh(x\v, x h) + Shh{xw, x h, 



C A {L^ H \x w ,x H ] 



6 ln(x 



/in / 



(2 S WH (x w , x H ) + S H h(x w , x h )) \ ■ 



(34) 



To obtain this, we have utilized Eq. (|T^) and the behaviour of Oll under Fierz-transformation. 
The finite quantities in front of the LO Inami-Lim functions in Eqs. (|33| ) and ( |34"D are the 
remnants of the functions Xa jfcO^o) °f Eq. (p4[). The coefficient Cll{^o) m Eq. (|29|) exhibits 
no dependence on the choice of the gauge and the external quark states. 

What remains to be done is the running of our effective theory down to a lower scale. In [|IJ it 
was shown that a good choice for this matching scale is /i = m w resulting in x M0 = 1 and that 
the physical observables do not depend on the choice of this scale. The scaling down from this 
matching scale to the mesonic scale is pursued by means of the renormalization group equation. 
The renormalization group equation for the Wilson coefficient Cll^) in the SM reads 



ClM = o 



(35) 



with the initial condition Cll(^o — m w) given in Eq. (|30|) . 

Expansion of the anomalous dimensions and the /3-function yields 



i{g) 



-A) 



Pi 



(4tt) 



(Ait) 2 
7 (0) + 7 (1) ^ 



(An) 



(An) 



+ ... 



(36) 
(37) 



The coefficient 7^ was calculated by Buras and Weisz in |13| . Note that 7 W is renormalization 
scheme dependent JL3, |19[ 2(J. The coefficients relevant in our further calculation are 



7 



7 



(i) 



1 



2n f ) 



3 (11 N ' 

A = ^Nl-^N c n f -2C F n f 
(o) = 6 ^-i 



N c 
N c -1 

2 N r 



57 19 , r 4 ' 
-21 H N c + - n f 

N r 3 3 i 



(38) 
(39) 
(40) 
(41) 



where N c is the number of colours and rif is the number of active flavours. 
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The solution to the renormalization group equation for the Wilson-coefficient reads 



ClM = exp 



S{mw) d ' 

m Pig*) 



C L L(m w ) 



(42) 



The evolution from the initial or matching scale mw to a lower one may be performed as in a 
massless theory because the top quark has been integrated out previously and appears in the 
matching condition only. At the NLO, Eq. fl42|) can be approximated by 



ClM ~ ex P 



exp 



g(mw) , 7 (0) ^/( 47r )2 + 7 (l) g^/^y 



dg 



C L L(m w ) 



7 



(0) 



+ 



7 



(i) 



«s(mw) 1 

«s(aO " 5 2 l/? ^ + A/(4vr) 2 ^ 2 ' /3 (4 7r) 2 + /3 1 ^ 



where we have used the Eqs. (^) and Q3"7|). After elementary integration, we obtain 



(43) 



1 + 



as(mwj ~ as(/i) / 7 (1) _ 7 (0) 
4vr bA) 2/3 2 



C L L(m w ) 



(44) 



where r] LO is the well known LO-scaling factor given in Eq. (|7]). 
In conclusion we find the following expression for H c r at NLO 



I Q2 

Hcs = 7^ m 2 w (V td V* b f 7] 2 {x w ,x H ) S mDM (x w ,x H )0 L L 

4 7T Z 



(45) 



with 



Oll 



1 | ots{m w ) ( D 2 }idm(xw,xhs z 



1 + 



4tt \S mDM (x w ,x H ) 



4 7T 







LL 



(46) 
(47) 



where 



7 



(i) 7 (o) 
2^1 



2/?o 



(48) 



Neither the factor r] 2 , nor the matrix element of Oll are dependent on the low energy scale 
ji « rriB, up to 0(a 2 s ). 
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4 Results 



We want to consider now the mass-splitting Am B between the electroweak mass eigenstates 
Bh and B^. The mass-splitting Amg is directly given by the S°S°-mixing amplitude 

Am B = —\(B°\H eS \W)\ 
m B 

G 2 

= T^ m w( v tdV t l) 2 S 2 HDM{xw,XH)V2( x w,XH) B B Jims ■ (49) 

We obtained the result in Eq. ( [49]) by using Eq. (g). B B is the renormalization prescription 
independent B-parameter. It is defined by U 



B B = B B (fi) a s (fi)' 6/23 



4tt 1 2 A, 2(3f\ 



(50) 



in 



The term proportional to as results from NLO-corrections. We use B B = 1.31 ± 0.03 [EI 
our numerical calculation. 

The current experimental mean value of the mass-splitting is given by Am B = (3.05 ±0.12) • 
10~ 13 GeV [p2] . We use the top-quark mass mj? olc = 175 ± 6 GeV |23[ which corresponds to 
m t (mw) = 176 GeV. Furthermore, we employ m w = 80.33 GeV, V t b = 0.9991 p4| in all plots. 



The meson decay constant is set to 175 ± 25 MeV |21|. 

From Fig. [5] we can deduce that the difference between the LO-calculation and the corre- 
sponding NLO-calculation for V t d is approximately 7 %, when mf ole is used in the LO expres- 
sion. Hence, the NLO-calculation enhances the result for the mass-splitting up to 15 %. This 
can also be concluded by a direct comparison of r] LO = as(M w ) 6 ^ 23 and r] 2 within the SM. We 
have found rj 2 = 0.4942 and r\ LO = 0.5751. 

If we consider an extended Higgs-sector within the framework of the 2HDM we have to 
decrease the related CKM-elements, to get an overlap between the allowed range for Am B and 
the Higgs-mass. We recognize for example in Fig. |6| that we can not find a physical Higgs- 
boson with a mass smaller than 1 TeV for V t d = 0.0086, assuming that the ratio of vacuum 
expectation values tan/? = 1 and f B \/B B = 0.2 GeV. If we decrease tan (3 we have to decrease 
Vtd even more. Fig. |5] is in full agreement with Fig. |6| in the SM limit m# — > oo. If we assume 
Vtd — 0.009 and if we furthermore take into account the errors for f B \/B B and m t the 2HDM 
can not be distinguished from the SM for Higgs-boson masses larger than 1 TeV. 

In Fig. 0, it is indicated that we have a relatively sensitive relation between the top-quark 
mass and the possible Higgs-boson mass. 

The comparison of NLO- and LO-calculation is shown in Fig. |^. We have plotted the mass- 
splitting over the Higgs-mass for two typical values of V t d and a ratio of the vacuum expectation 
values tan/5 = 1.25. The difference between NLO- and LO-calculation approximately amounts 
to 18 %. 
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5 Conclusions 



We have calculated the 5°-B°-mixing within Next-to-Leading-Order with the inclusion of two 
charged Higgs-bosons. The NLO calculation leads to an effect of approximately 7 % for Vu in 
comparison to the LO calculation within the SM. This is indicated in Fig. [|. 



We have verified that the scheme developed by Buras, Jamin, Weisz [[y, [L3| for a proper 
separation of long and short distance effects in QCD remains valid in our 2HDM calculation. 
This can be considered as a good cross-check of the calculations presented here. The inclusion 
of NLO-corrections lowers the mass-splitting in the SM as well as in the 2HDM. We find, for 
example, a difference between the pure LO- and the NLO-calculations of approximately 18 % 
for rrtu — 200 GeV and 17.5 % for = 400 GeV. Hence the NLO-contributions plays an 
important role in the 2HDM as they correct the LO result by about 18 %. 

The inclusion of the two Higgs-bosons leads to an increase of the calculated amplitude, which 
depends on the mass and the vacuum expectation values of two Higgs doublets. In all practical 
calculations, the lower bound of the Higgs-boson mass was assumed to be 100 GeV and the 
upper bound was 1 TeV. 

To obtain a mass-splitting within the experimental allowed region, it is necessary to decrease 
the CKM-matrix elements for small values of tan/3 and m#. Our calculations are not valid for 
higher values of tan/3 (region near tan/3 = 40), because it is then necessary to introduce new 
operators. 

In the limit of very heavy Higgs-bosons we verify the well-known results for the Standard 
Model. 
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Appendix A 

A.l Special functions 

After carrying out the second integration, the dilogarithm or Spence function appears in our 
results. It is defined by 

Li 2 (x) = - rd*k(lz*)=f;^ ; | X |<1. (51) 
Jo t ^ n 2 

The following useful relations hold 

Li 2 (l-x) + Li 2 (l-l/x) = -^ln 2 (x) , (52) 
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7T 2 

Li 2 (x) + Li 2 (1 - x) = — - ln(x) ln(l - x) . (53) 

A. 2 Inami— Lim functions and NLO in the SM 

We list now the Inami-Lim functions and the well-known functions L^ l,SM \xw)- 

(I 9 3 3x 2 w \n(x w )\ 

o ww (x w ) — x w [-7 + ~n-\ V ~~ oTi \2 _ oTi \3 ' ^ 

\4: 4(1 — xw) 2(1 — x w ) 2(1 — xw) J 



2 SVh^vk, — 



1 xh xw ( 1 + 2xh ln(x_ff) 

t;an%9) 1^(1 - x H f + (1 - x H ) 3 

1 a;^ x w / (2x w — 8x#) hi(x#) 6x H /ln(x w ) 



Shh(xh) = —47^ — r~ ( - — + — ) , (55) 



tan 2 (/5) 4 \(1 - x H ) 2 (x H - x w ) (1 - x w ) 2 (x H - x w ) 
8 — 2xw \ 

(1 - X H )(l - X W ) J 



(56) 



The function Swh {xw, x h) is obtained by including the charged would-be Goldstone bosons or 
in a more direct way, by choosing the unitary gauge for the W^-boson right from the beginning. 
The functions L^' SM ^ can be decomposed as 

L^ SM \x w ) = WW®(x w ) + 2W$W(x w ) + (57) 
where 

WW^ixw) = WW t { p(x w )-2Ww£\xw) + WWV(x w ) , (58) 



\ (4 x w + 38 x 2 w + 6 x* w ) \n(x w ) , (12 x w + 48 x^ + 12 x^) Li 2 (l - l/x w ) 
W vv u (Xw) — : —4 1 : 7-4 



(x w - 1) (x w - 1) 

(24 x w + 48 xfv) Li 2 (l - x w ) 3 + 28x^ + 17x2^ 



(x w - l) 4 (x w ~ I) 3 



2wwW( Xw -) = 2(3 + 13 x w ) 2x w (5 + 11 x w ) \n(x w ) 



(59) 



(x w - l) 2 (x w - l) 3 

12x^ (1 + 3xw) Li 2 (l — l/xw) 2Ax w (1 + x w )\j\ 2 (1 — x w ) 



(x w - l) 3 (x w - 1) J 



(60) 



WW<£(x w ) = 3 , (61) 
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x w (7 + 52 x w - 11 x w ) 3 x 3 ^ (4 + 5 x w - x w ) \yy(x w ) 



A(x w -lf 2(x w -iy 
3x w (3 + Ax w - x^)Li 2 (l - l/x w ) 18 xfy Li 2 (l - x w ) 

4x^ (11 + 13 xw) 2x w (5 + x^)(l — 9xw) ^(xw) 
(x w - I) 3 + - I) 4 

24x^(1+ 4 x^/ + x^)Li 2 (l - 1/%) 48 ^ (1 + 2x w )U 2 (l - x w ) 



(xw - I) 4 (xw - I) 4 

(63) 



WW t ( P (x w )-2 WWg> (x w ) + WW™ M ; ( 64 ) 
2 x w (A — 3 x w ) \n(x w ) (12 x^ — 12 x w — 8 xf^) Li 2 (l — l/xw) 

(65) 



(8 - 12 xvk + 12 x w ) Li 2 (l - x w ) (23 - x w ) 



(xw - I) 4 (xw ~ I) 2 



2 (2 — xvk) tt 2 (8 — 5 x^) ln(xw) 



3a?w (x w -l) 
(& x w + A x w ) Li 2 (l - l/x w ) (8 + 12 x w - 6 x w ) Li 2 (l - x w ) 



+ 



x w (xw — I) 2 x w (xw — I) 2 
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(x w - 1) 
4 



(66) 



-23 + - tt 2 , (67) 



11 x^, (1 + xw) x 3 ^ (4 - 3 x w ) hi(x w ) 



A(xw-l) 2 2(x w -l) 3 
x w (3 - 3x w + 2x^)Li 2 (l - l/x w ) x w (2 + 3 x w - 3 x^)Li 2 (l - x w ) 
(xw-lY + (xw-l) 4 

(61 

12 



{8) 30 x 2 w \2x z w ln(xty) 12x^Li 2 (l - l/x w ) 

2W * M = w^w + ixw-if w^w — 

I2x 2 w {2 - x 2 w )U 2 {l-x w ) 

W^W ' (69) 

A. 3 NLO in the extended Higgs sector 

The L^' H ^ can be decomposed in the following way 

i( *' ff w,*») + tJ&sj M<,, (— ») + tJm hh " Hxh) ■ 

(70) 

We list now the newly obtained results for the explicit QCD corrections in the extended Higgs 
sector: 



HHW( Xi ) = ^^ (1 \x H ) + 6(Hx H )-Hx w )) xr£- S HB (xi) , (71) 

%H i=H,W aXi 



(1) / 2a:^(13 + 3a: g )ln(a; g ) 2 x H (9 + 7 x H + 7 x w -23x w x H ) 

[Xi) ~ XW \{x H -lf{x H -x w ) (x w -iy( XH -iy 

2x^(18 — 6x H — AAx w + 13 x H x w + 9x H x^) \n(x w ) 

(x H - l) 2 (x w - l) 3 (x H - x w ) 
2x H x w {§ — 27 x w + Gxlv + Qx H Xw) ^ n ( x w) 2Ax 2 H \n(x H ) \n(x w ) 



(x H - l) 2 (x w ~ l) 3 (x H - x w ) (x H - l) 3 (x H - x w ) 

24x^Li 2 (l - l/x H ) 2Ax H x w (1 + x w ) Li 2 (l - l/x w ) 



(x H - l) 2 (x H - x w ) (x w - l) 3 (x H - x w ) 

ASx w x H Li 2 (l -x w ) \ 
(x w - lf(x H -x w ) ) 



9 ( b TT(i)( \ 2 / x H (31 -15x H - 15 x w - x H x w ) x h (7 + 21 x H - 12 x 2 H ) ln(x H ) 

[Xi ) x w y 2{ ^_ i)2 _ 1)2 2 { ^ _ 1)3 _ ^ 

xh (7 — 9x w + 36 

X\y 

18 x^) \n(x w ) 
2(x H - l) 2 (x H - x w )(x w - l) 3 
_^x 2 H (% - 36x w + Qxly + 3x 3 v )\n(x w ) 
(x H - l) 2 (x H - x w )(x w - l) 3 
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x H (11 — 45 x w + 18 x^) ln(xv^) ^ 6x# ln(x#) ln(x^) 

2(x# - l) 2 (x# - Xn/)(x W - l) 3 - l) 3 (x H - x w ) 
Gx H (l + x H - x 2 H )U 2 (l - IJxh) 

(x h - l) 2 (x H - x w ) 
6x H {l + 2 

X\y ~ ^v^)Li2(l — l/#w) 12 X# Ll2(l — %) 

(% - xhO(x^ - I) 3 (xh - %)(% - I) 3 



^$$( 8 )(^) + 6 (ln(^) - ln(x w )) S HH (x t ) , 



\(x# - %)(% - l) 2 
6x^(5x w - x H + 3x^x^^2(1 - l/x w ) 

(x H - 1) 2 (xh - %)(% - l) 2 x w 
6x H (2x w - lOx H x w + x H x w )Li 2 (l - 1/W) 

(x H - l) 2 (x H - x w )(x w - I) 2 
6x 2 H (5x w -x H - 8x w + 2x H -x 2 y )Li 2 (l - x H ) 

(x h - l) 2 (x H - x w ){xw ~ ^) 2 x w 
Q(xly - x H x w + 2x 2 H Xw)Li 2 (l - x H ) 

(x H - l) 2 (x H - x w )(x w - I) 2 
6x 2 H (-x H + 5xjy)Li 2 (l - l/x H ) 

(x H - l) 2 {x H - x w )x w 
Qx 2 H (5x w - x H - 8x w + 2x jH -x 2 y )Li 2 (l - x H /x w ) 

(xh — l) 2 (x// — x w )(xw — l) 2 xw 
6(x w - x H x w + 2x^x^)112(1 - xh/xw) 

(x h - l) 2 (x H - x w )(x w - l) 2 
6x g (l -x H - ln(xg)) 6x H (2x w - 1) ln(arw) 

(x H - l) 2 (x w - 1) (xh - l)(xw - I) 2 

Qx 2 H (5xw — xh — 8x w ) ln(xn) ln(x^) 

(x H - 1) 2 (xh - x w )(x w ~ ^) 2 xw 
l2x 2 H (x H x w + Xw) \n(x H ) \n(x w ) \ 
(x H - l) 2 (x# - x w )(x w - I) 2 / 



\ 2(x H - 1) (x w - 1) x w 
(24 - 7 x H x w + 2 x 2 H x w + 2 x^ + x H x 2 w ) ln(x g ) 

2 (% - l) 2 (x H - X^y) (x W - 1) X W 
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xh (7 - 7 x H + 4 x w - 6 x 2 ^) ln^iy) 
2 (% - 1) (x H - £w) (^w - I) 2 
| (x 2 ^ + x 2 w - 3 x|- x 2 ^) \n(x w ) 

(X H ~ 1) (^ff - X W ) (X W - l) 2 X W 

x\ (4 - 6 x w + 3i H %) ln(xH) hi(xvy) 
(x// - l) 2 (x# - x w ) (x w - l) 2 x w 

- 3x^ + 6x 3 n, - 3 Xw) \n(x H ) hi(x w ) 
(x H - l) 2 (x H - x w ) (x w - l) 2 x 2 w 

XH (3 X w + 2 X g Xiy (2 + - X^ (1 + 2 Li 2 (l - l/x H ) 

(x H - l) 2 (s H - x w ) x 2 w 
(4 x H x w - 6 x 2 ^ x w + 3 x 2 H x 2 w - x 2 y )Li 2 (l - x H ) 

(x H - l) 2 (x H - x w ) (xw - I) 2 x H 
(4x 2 H x w - 6 x 2 H x^y - x 3 H + 3 x|- x 2 y )Li 2 (l - x H ) 
(x H - l) 2 (x# - x w ) (x w - l) 2 x 2 w 
^ 2 x 2 H (6 - x w - 3 x H + x w x H ) Li 2 (l - l/x w ) 
(x H - l) 2 (x H - x w ) (x w - I) 2 
xh (3 xly + 4 % x w - x 2 ^) Li 2 (l - l/x w ) 
(x H - l) 2 (x H - x w ) (x w - l) 2 xly 
^ (4x H x w - 6 x 2 H x w + 3igXy - x 2 v )Li 2 {l - x H /x w ) 
(x H - l) 2 (x H - x w ) (x w - I) 2 x H 
xM 4 % - 6 x w - x H + 3 x H x w )Li 2 (l - x H /x w ) 
(x H - l) 2 (x H - x w ) {xw - I) 2 x w 
6 x g Li 2 (l - x^) \ ^ 76 ^ 
(x H - x w ) (x w ~ I) 2 / 

We have used the capital W, $ and H to label the W-boson, the charged would-be Goldstone 
boson and the physical Higgs-boson, respectively, i.e. denotes a diagram with one internal 
W^-boson and one internal unphysical scalar. 
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Figure 1: Box diagrams for the i3°i3 -mixing in the framework of 2HDM. We have also taken into account 
crossed diagrams, which are related to the original ones by a Fierz-transformation. 







Figure 2: QCD-corrections to the effective four quark interaction. They are needed in the calculation of the 
anomalous dimension at LO and for the correct separation of long and short distance contributions. 
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Figure 3: Diagrams for the NLO-calculation in the SM, which are convergent in the limit of vanishing external 
masses. Additional diagrams with one Higgs-boson and two Higgs-bosons have to be calculated in the 2HDM. 
Furthermore, we have to consider crossed diagrams. 
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Figure 4: Diagrams for the NLO-calculation in the SM, in which infrared divergences appear when one sets 
external masses to zero. The last diagram is the so-called double-penguin, which does not contribute to the 
mixing at the leading order in (mj /m\y) 2 ■ Additional diagrams with one Higgs-boson and two Higgs-bosons 
have to be calculated in the 2HDM. Furthermore, we have to consider crossed diagrams. 
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Figure 5: The CKM-element Vtd is plotted versus the B meson-decay constant times the square root of the 
B-parametcr. The shaded areas represent the allowed values when errors of Am^ and rat are included. There 
is a small overlap between the LO- and the NLO region shown in the picture. 
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Figure 6: The mass splitting within the 2HDM including NLO-corrections for different values of Vtd and for 
tan/3 = 1. The top-quark mass is set to ra\ aXc = 175 GeV. The shaded strip is the experimentally allowed 

1 /2 

region. The factor fsBg is fixed to 0.2 GeV in this figure. The limit ran — > oo yields the Standard Model 
results. 
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Fi gure 7: The mass splitting within the 2HDM for different values of the top-quark pole mass. In this figure, 
V td = 0.009, tan/3 = 1, fsB^ 2 = 0.2 GeV. 
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Figure 8: A comparison between LO calculations and NLO-calculations within the framework of the 2HDM. 
Here, tan/3 = 1.25, f B B^ 2 = 0.2 GeV and m? t olc = 175 GeV. 
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